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, , We discuss the cyclic Wilson loop, i.e. a rectangular Wilson loop that spans the entire compacti- 

fled time axis in the imaginary time formalism. The result of a perturbative calculation at 
is given, with the main focus on the ultraviolet divergences of this operator. Based on a general 
analysis of divergences in loop diagrams, we find that, unlike usual Wilson loops, the cyclic loop 
does not have cusp divergences but intersection divergences, where the intersections arise due to 
the periodic boundary conditions. As a result, the cyclic Wilson loop itself is not multiplicatively 
^ renormalizable, but the difference between it and the correlator of two Polyakov loops is. 
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1. Introduction 

Wilson loops are very useful operators to study heavy quarkonium systems [1, 2, 3]. In the 
infinite mass limit the propagator of a heavy quark can be expressed through a straight Wilson line 
in the time direction. Wilson lines also allow for the definition of gauge invariant states of two 
spatially separated quarks. Therefore the time evolution of a gauge invariant quark-antiquark state 
can be described by the expectation value of the trace of a rectangular Wilson loop. In the infinite 
time limit one can project out the lowest energy eigenvalue [4, 5], and it's possible to account 
for finite mass corrections through insertions of gluonic operators into the Wilson loop [2]. At 
finite temperature there is a similar operator called the Polyakov loop, which consists of a straight 
Wilson line along the whole of the compactified time axis in the imaginary time formalism. The 
expectation value of the correlator of two traced Polyakov loops is related to the free energy of 
a system of an infinitely heavy quark and antiquark, and has a known spectral decomposition in 
terms of the energy eigenvalues of the static QCD Hamiltonian [6, 7]. A Wilson loop, called the 
cyclic Wilson loop, with two Polyakov loops as its time-like sides, is expected to have similar 
interpretations. 

This paper is concerned with the properties of the cyclic Wilson loop, in particular its diver- 
gence structure. Typically, rectangular Wilson loops have ultraviolet divergences in addition to the 
usual ones related to charge renormalization. These additional divergences are called cusp diver- 
gences, because they are related to the four right-angled corners of the loop's contour. They are 
known to factorize and can be removed by a multiplicative constant [8, 9]. However, it was found 
that the additional divergences of the cyclic Wilson loop are different and cannot be removed by a 
multiplicative constant [10]. We present a solution to this problem by investigating the origins of 
UV divergences in terms of loop diagrams, and demonstrate that the difference between the cyclic 
Wilson loop and the correlator of two Polyakov loops is indeed multiplicatively renormalizable. 

The paper is organized as follows. Section 2 will detail some properties of the cyclic Wilson 
loop, which will be instrumental in the following discussion, and give a perturbative calculation 
at &(gCs). Section 3 will then discuss origins and different types of UV divergences in loop func- 
tions and apply this analysis to the case of the cyclic Wilson loop. Section 4 will then test the 
general statement of section 3 at ff(a^), while section 5 contains our conclusions. All these results 
are explained in more detail in our paper [11]. 

2. Properties of the cyclic Wilson loop 

The Wilson line operator, for a line going from y to x along the contour C, is given by 



where & stands for path ordering of the colour matrices in (z) along the contour C. Since we will 
use the imaginary time formalism of finite temperature field theory, we work in Euclidean space 
with periodic boundary conditions in the time direction. Under a gauge transformation V G SU (N) 
Wilson lines transform as V (x)W (x,y;C)V^ (y), so that the trace of a Wilson line with a closed 
contour, i.e. a Wilson loop, is a gauge invariant operator. 




(2.1) 
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A Polyakov loop operator is given by a straight Wilson line, which runs once around the 
compactified time axis: 

r r l/T 

ig / dTA (T,r) . 
Jo 



P(r,T) = ^exp 



(2.2) 



Because of the periodic boundary conditions, the contour of a Polyakov loop is closed and the trace 
over a Polyakov loop is gauge invariant. The correlator of two Polyakov loops, hereafter simply 
called the Polyakov loop correlator, is then given by 

P c (r,T) = (fr[P\0,T)}fr[P(r,T)}) , (2.3) 

where Tr denotes the trace over colour matrices divided by the number of colours N. A single 
Polyakov loop operator does not depend on the position vector r because of translation invariance, 
but the correlator of two Polyakov loops does depend on their distance r = |r|. P c is related to the 
free energy of an infinitely heavy quark and antiquark in a thermal medium. 

The cyclic Wilson loop is defined similarly to the Polyakov loop correlator, but instead of 
having two separate traces over the Polyakov loop operators, their endpoints are connected by two 
spatial Wilson lines to obtain a gauge invariant object. The cyclic Wilson loop therefore consists 
of four straight Wilson lines, the two corresponding to Polyakov loops will be called quark lines, 
the two spatial lines will be called strings. 

W c {r,T) = (fr[W(r,0)P\0,T)W(0,r)P(r,T)]) . (2.4) 

In the arguments of the string operators only the spatial components are displayed, the temporal 
components are always or 1/T, which is equivalent. Because Wilson lines are unitary operators, 
and because exchanging start and end point of a Wilson line corresponds to taking the hermitian 
conjugate, the two string operators are inverse to each other. This leads to many cancellations of 
Feynman diagrams in a perturbative expansion of the cyclic Wilson loop. Since this is a direct 
consequence of the periodic boundary conditions - in a usual rectangular Wilson loop the strings 
are separated in the time component and therefore not inverse to each other - we will call this effect 
cyclicity cancellation. 

A perturbative expansion of the cyclic Wilson loop gives Feynman diagrams, which consist of 
the rectangular contour specified in (2.4) and gluons attached to points on this contour. These end- 
points of gluon propagators each give a factor of ig and are integrated along the contour according 
to the path ordering prescription (cf. eq. (2.1)). We will refer to them as line vertices, as opposed to 
internal vertices connecting gluons to other gluon, ghost or quark propagators (cf. ref. [12]). Now, 
cyclicity cancellation means that all diagrams cancel which have at least one line vertex on a string 
and no line vertex on at least one quark line. 

Another useful feature of the cyclic Wilson loop (or in general any closed Wilson loop) is 
expressed in the exponentiation theorem [13, 14], which says that the perturbative expansion in 
Feynman diagrams can be resummed and exponentiated: 



W c = ^C(y)W(r)=exp 



I C( 7 )W(7) 
ye2PI 



(2.5) 
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where y denotes the Feynman diagrams, C(y) represents the colour factor of a diagram, which is 
given by the normalized trace over the colour matrices and contractions with structure constants 
from this diagram, and W(y) the value of a diagram without its colour factor, i.e. the result of 
applying the Feynman rules to a diagram and performing the corresponding integrations. Eq. (2.5) 
states that the logarithm of the cyclic Wilson loop can be expressed as a sum of the same Feynman 
diagrams, but only the subset of two-particle irreducible (2PI) diagrams contributes. A diagram is 
reducible, if the contour can be cut in two places such that the two pieces of the contour are not 
connected to each other by gluon, ghost or light quark propagators (excluding the case where there 
is no gluon attached to one of the pieces at all). If this is not possible, then it is a 2PI diagram. So 
the logarithm of W c is given by a series over only 2PI diagrams, however their colour factors have 
to be replaced by the colour-connected coefficients of each diagram. These are denoted by C in 
eq. (2.5). We refer to refs. [13, 14, 11] for an explicit definition of these coefficients. 

The advantage of the exponentiation theorem is, that it is in general easier to calculate colour- 
connected coefficients instead of calculating reducible diagrams. So in the following we will need 
to consider only 2PI diagrams, and furthermore only those 2PI diagrams which do not cancel 
through cyclicity. For calculations, which have been regularized in dimensional regularization with 
d = 3 — 2s, Coulomb gauge turns out to be very advantageous, because in this gauge several of the 
remaining diagrams simply give zero. The relevant diagrams at &(ccf) are given in figure 1. After 
removing the UV divergence from the gluon self energy through charge renormalization (in MS, 
with 1/e = 1/e — Ye +ln47i) the result reads 



InWr 




31 20 \ /ll 4 \ 
yC A - — T F n f \ + I —C A - -T F n f \ Inplr 

f d^k (e' r ^ l) 1 

+ i%a - r I (2,1)3 f (r 11 ™ (°' k >) I + c fCa^ 

+ Cj £r (i + 1 +2ft -^ + lnf V + 2£ + . (2.6) 

The square brackets in the first line come from the vacuum part of the gluon self energy (e.g. [15]), 
the second line starts with the contribution from the thermal part [16, 17], whose integral is not 
performed analytically but is UV and IR finite. The other term in the second line comes from 
the diagrams with a 3-gluon vertex. The last line gives the contribution from the last diagram 
in figure 1, where the series expansion of the thermal part is valid for rT < 1. This is the only 
divergent contribution at this order in a s . The integration for the gluon propagator connecting the 
two quark lines in this diagram has been carried out in 3 dimensions. If this integral had been 
computed in d dimensions, then its 0(e) part would give a finite contribution when multiplied by 
the 1 /e term in eq. (2.6). But after the divergence is subtracted as shown in the following section, 
this finite term will also be removed, so it is not necessary to display it here. 

We have implicitly assumed the hierarchy of scales 1/r ^> T S> mo 3> ct s /r here, and give 
only the contribution from the highest scale 1/r. This discussion deals with UV divergences, so 
the low-energy behaviour of finite temperature QCD is not relevant. The general discussion of UV 
divergences in the following section, however, will not depend on this hierarchy of scales. The 
result (2.6) agrees with the appropriate limit of the calculation found in [10]. 
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Figure 2: Illustrations for the different types of loop function divergences. 



This result shows the peculiar difference between W c and a usual rectangular Wilson loop with 
cusps regarding UV divergences. There is no cusp divergence at ${a s ) (it has cancelled because 
of cyclicity), instead the divergence at 0{ccj) is not a constant term, but multiplies a function of r 
and T. It can therefore not be removed through a multiplicative constant to W c , which would be an 
additive constant to eq. (2.6) because of the logarithm. The reason for this will become apparent in 
the next section. 

3. UV divergences in loop functions 

Here we will first discuss the origin of UV divergences in loop functions, i.e. expectation val- 
ues of traced Wilson loops, and then apply the findings to the cyclic Wilson loop. This first part 
is based on the analysis found in [8, 9]. In position space, UV divergences arise when vertices 
are contracted at one point. In the case of internal vertices this leads to the usual UV divergences 
removed by charge renormalization. For loop functions it is also possible to contract line vertices, 
giving rise to a new sort of UV divergences. Depending on where and how many of these line ver- 
tices are contracted, we can further distinguish different types of UV divergences. The superficial 
degree of divergence for such a contraction is given by 

<W sm ooth = 1 ~N ex ^singular = ~N ex . (3.1) 

iVex denotes the number of external lines, i.e. propagators leading to vertices which are not at the 
contraction point. The subscript refers to the point of contraction, where the contour can be smooth, 
have an angle or intersect itself. The latter two cases are called singular points. Then there are three 
possibilities: (a) line vertices are contracted at a smooth point without external lines, which gives 
a linear divergence, (b) line vertices are contracted at a smooth point with one external line, which 
gives a logarithmic divergence, (c) line vertices are contracted at a singular point without external 
lines, which also gives a logarithmic divergence. 
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Linear divergences, like those in (a), are automatically removed in dimensional regularization. 
They are therefore not relevant to this discussion. In other regularizations they may be absorbed in 
a "loop mass" term, which is proportional to the length of the contour. 

The logarithmic divergences (b) will be called line vertex divergences, because they give in 
effect a correction to the line vertex of the external propagator and will be removed by a line 
vertex counterterm. This counterterm comes in through the factor ig at each line vertex and is 
automatically included in the charge renormalization procedure. This was shown in [8]. In the 
result (2.6) presented above charge renormalization has already been carried out, so line vertex 
divergences cannot be responsible for the divergence found. 

The logarithmic divergences (c) will be called cusp or intersection divergences, depending on 
the nature of the singular point. It was shown in [9] that cusp divergences factorize out of loop func- 
tions and can therefore be removed through a multiplicative constant. This constant only depends 
on the cusp angle. Intersection divergences were also considered in [9], where the authors showed 
that loop functions with intersections mix under renormalization with other loop functions which 
have identical contours apart from a different path ordering prescription at the intersection. So for 
each intersection point there is a renormalization matrix mixing these associated loop functions. 
The renormalization matrices depend only on the angles at the respective intersection. 

All four cases are illustrated in fig. 2. The dashed boxes mean that all vertices inside the box 
are contracted at a point on the contour. If the box contains a singular point, then this is implied 
to be the contraction point. The first two diagrams are examples for linear divergences. The next 
two lead to a line vertex divergence which is removed through the counterterm in the following 
diagram. The last two diagrams in the first row are examples of cusp divergences. The dependence 
of the divergence on the cusp angle is introduced through diagrams like these, similar diagrams with 
all line vertices before or after the cusp also contribute to the cusp divergence but do not depend 
on the cusp angle. The second line shows the contours of two associated loop functions with a 
self-intersection, which have to be renormalized together. The contours are identical everywhere 
except for the intersection. If we denote the left and right loops as Wi and Wr respectively, then 
the first loop function is given by (Tr [WlWr]) and the second by (TrfW^jTrfW/j]). The last two 
diagrams are examples for intersection divergences. The first one of these appears in both loop 
functions and is identical to the contribution of a corresponding cusp divergence. The second 
example only contributes in the case of the first contour. It is diagrams like this which introduce 
the dependence on all intersection angles and necessitate the mixing of associated loop functions 
in order to renormalize them. 

Turning back to the cyclic Wilson loop, due to the periodic boundary conditions the singular 
points are no longer cusps but intersections, because the two strings run along the same points 
in Euclidean space-time. In principle all points on the strings are intersections, but it can be 
shown that intersections where all angles are either or n are not divergent, so only the two 
string endpoints need to be considered as intersections. The associated loop functions are given 
by (Jr[W(r,0)P t (0,r)W(0,r)]fr[P(r,r)]^) and (Tr [P f (0,T)) fr [W(0,r)P(r,T)W(r,0)]j for 
choosing a different path ordering at the upper or lower intersection respectively, and for choosing 
a different path ordering at both intersections by (ft [P^(0,T)] Tr[W(r,0)W(0,r)]Tr[P(r,7 , )]y 
In all three cases the string operators cancel, so every associated loop function corresponds to the 
Polyakov loop correlator. 
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So cyclic Wilson loop and Polyakov loop correlator are renormalized together by a 2 x 2 
matrix. Observing that the Polyakov loop correlator is already finite after charge renormalization 
and that the angles at both intersection points are identical is enough to derive the following form 
of this renormalization matrix: 



w c {Ry 




(3.2) 



By rearranging this equation we see that W c — P c is multiplicatively renormalizable, while W c 
alone is not. Using eq. (2.6) and the fact that P c = l + &(g 3 ) (e.g. [18, 19]), we can determine the 
renormalization constant at first order to be Z = 1 — S^ 1 + €?(aj). The renormalized Wilson loop 
is then given by eq. (2.6) with the only modification that the l/e term is absent. 

4. Test of the renormalization formula at &(af) 

In order to provide a test for this renormalization formula (3.2), we will discuss the intersection 
divergences arising at &{a 3 ). Since the first order term of Z = 1 +LnZi° ; " has already been de- 
termined by the 0(ocj) calculation, Z\ = — 1|, all contributions to the renormalized cyclic Wilson 
loop at &(ots) involving Z\ must cancel in order for eq. (3.2) to be valid. 

We will start by considering which diagrams can in principle contribute to the intersection 
divergence. The restriction to 2PI diagrams which do not cancel through cyclicity leaves only two 
possibilities for the remaining diagrams: either they have all their line vertices on the same quark 
line and none on the strings, or they have some line vertices on both quark lines. In the first case, 
these diagrams contribute equally to the Polyakov loop, and since the Polyakov loop is finite after 
charge renormalization, they cannot contribute to the intersection divergence of the cyclic Wilson 
loop. So divergent diagrams need to have at least one gluon connecting one quark line to the 
other. As explained in the previous section, intersection divergences arise when line vertices are 
contracted at the intersection without external lines. The gluon(s) connecting the quark lines cannot 
be contracted at any of the intersections, so there can be no connections between these gluons and 
the line vertices contracted at the intersection, otherwise there would be external lines. 

So we can think of the intersection divergent diagrams as having two components. One con- 
sists of the gluons connecting the quark lines, which by themselves would be a finite diagram. We 
will call this the base diagram. The other component consists of additional gluons added to the 
base, in such a way that they can be contracted at the intersection without external lines. These 
additional gluons need to have line vertices to the left and the right of the ones of the base diagram, 
otherwise this diagram would be reducible. 

So in order to find all intersection divergent diagrams at @{ol 3 ) we just need to take all lower 
order bases and add gluons that can be contracted at the intersection. All relevant bases are shown 
in fig. 3. Again we will use Coulomb gauge for simplicity, otherwise there would be more diagrams 
to be considered. A treatment in more general gauges can be found in ref. [11]. 

Fig. 4 shows intersection divergent diagrams with the first diagram of fig. 3 as base. Not 
all divergent diagrams with this base are displayed, but the remaining ones can be obtained from 
these by taking the minor image or moving some line vertices from the strings to a quark line (but 
not beyond the line vertices of the base). The sum in fig. 4 is supposed to include all of them. 
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Figure 3: A list of all relevant bases up to 0{(X% ) in Coulomb gauge. The last diagram is not 2PI, but after 
adding a gluon to make it intersection divergent it will be 2PI. 



Z\a s x 




= -Z 2 aix 



Figure 4: A sample of diagrams that determine the value of Z2. The dots stand for remaining finite terms. 
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Figure 5: Intersection divergent diagrams which are renormalized similarly to the (?(a,f) divergence. Again, 
the dots stand for remaining finite terms. 
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Figure 6: The cancellation of the remaining divergences in W c — P c at G{(X^) involves the Polyakov loop 
correlator. 



The divergences are proportional to a s 2 , so they have to be cancelled by Z2. By calculating these 
diagrams, which we will not do here, one can determine the value of Z2. 

The other bases in fig. 3 all need one additional gluon connecting the strings to become inter- 
section divergent &{o^) diagrams. Except for the last one, the cancellation of their divergences is 
very similar to the case at ff(aj). The added gluon gives the same divergent factor multiplying the 
base diagram and the divergence gets cancelled by Z\ times the base. This is illustrated in fig. 5. 

The cancellation of divergences for diagrams with the last base in fig. 3 is more complicated. 
An illustration is given in fig. 6. Since the sum over the diagrams shown in 1 gives the logarithm 
of the cyclic Wilson loop, we need to re-expand the exponential when we calculate W c — P c . This is 
where the product of diagrams on the left-hand side and the square of the first diagram on the right- 
hand side come from. The last diagram in fig. 6 comes from the Polyakov loop correlator, whose 
value at ff(aj) is given by P c = 1 - (Cf - \C F C A )^f2 + 28P + 0(g 5 ) (e.g. [18, 19]). The part 5P 
comes from diagrams common to both W c and P c , like the third and fourth diagram of fig. 1, and 
therefore cancels in W c — P c . The relation shown in fig. 6 in particular confirms the renormalization 

(R) 

formula W c y ' - P c = Z(W C -P c ), since here also the first non-trivial expansion term of P c plays 
an important role in the cancellation of divergences. This exhausts all possibilities for intersection 
divergent diagrams at ^(ot s 3 ) in Coulomb gauge. 
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5. Conclusions 

The cyclic Wilson loop has been introduced and calculated in perturbation theory at G(a^). 
Using a general study of the divergences of loop functions we have analyzed the divergences of 
the cyclic Wilson loop and found that they do not originate from cusps but from intersections. 
Accordingly the cyclic Wilson loop cannot be renormalized by itself, one rather has to mix it with 
the Polyakov loop correlator to obtain a finite quantity. As a result of this approach we find that 
W c — P c is in fact a multiplicatively renormalizable operator. This relation is confirmed by a study 
of the intersection divergent diagrams at ^(a 4 3 ). We have restricted ourselves to Coulomb gauge 
here, a generalization to other gauges as well as other extensions of this work can be found in [1 1]. 
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